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Abstract 


This  investigation  is  concerned  with  the  dynamic  simulation  of  com- 
plex structures  consisting  of  an  assemblage  of  substructures.  Each  of 
the  substructures  possesses  a large  number  of  degrees  of  freedom.  In 
fact,  in  theory  a continuous  substructure  has  an  infinite  number  of  de- 
grees of  freedom,  but  for  practical  reasons  each  of  the  substructures 
must  be  simulated  by  a limited  number  of  degrees  of  freedom.  The  com- 
ponent mode  synthesis  method  advocates  representing  the  substructure 
motion  by  a given  number  of  substructure  modes.  In  an  earlier  paper, 
the  authors  of  this  report  have  argued  that  component  modes  are  actually 
not  necessary  and  admissible  functions  suffice.  Equally  important  is 
the  treatment  of  discrete  substructures  with  a large  number  of  degrees 
of  freedom.  To  simulate  the  motion  of  such  substructures,  these  authors 
are  advancing  the  concept  of  "admissible  vectors".  The  first  phase  of 
the  research  has  concentrated  on  the  development  of  this  concept  and 
placing  it  on  a sound  mathematical  foundation  as  well  as  developing 
methods  for  producing  such  vectors.  To  introduce  the  above  ideas,  a 
discrete  model  consisting  of  a rotating  lumped-parameter  cantilever  beam 
has  been  used.  Such  a model  is  of  special  Interest,  as  it  can  be  re- 
garded as  simulating  a helicopter  blade. 
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1.  Introduction 

This  investigation  is  concerned  with  the  dynamic  simulation  of  com- 
plex structures  consisting  of  an  assemblage  of  substructures.  Each  of 
the  substructures  possesses  a large  number  of  degrees  of  freedom. 

In  fact,  in  theory  a continuous  substructure  has  an  infinite  number  of 
degrees  of  freedom,  but  for  practical  reasons  each  of  the  substructures 
must  be  simulated  by  a limited  number  of  degrees  of  freedom.  The  com- 
ponent mode  synthesis  method  (Ref.  1)  advocates  representing  the  sub- 
structure motion  by  a given  number  of  substructure  modes.  This  seems  to 
point  to  the  necessity  of  producing  a suitable  set  of  modes,  perhaps  by 
solving  the  eigenvalue  problem  associated  with  the  substructure.  Quite 
often,  however,  this  is  no  easy  matter  and  in  some  cases  the  eigenvalue 
problem  cannot  even  be  defined  properly. 

In  an  earlier  paper,  the  authors  of  this  report  have  argued  that  com- 
ponent modes  are  actually  not  necessary  to  represent  a substructure  mo- 
tion. Indeed,  if  the  substructure  is  represented  by  a distributed-para- 
meter mathematical  model,  then  admissible  functions  suffice  (Refs.  2-4). 
Equally  important  is  the  treatment  of  discrete  substructures  with  a large 
number  of  degrees  of  freedom.  To  simulate  the  motion  of  such  substructures, 
these  authors  are  advancing  the  concept  of  "admissible  vectors". 

The  first  phase  of  the' research  work  has  concentrated  on  the  develop- 
ment of  the  concept  of  admissible  vectors  for  discrete  systems,  placing 
the  concept  on  a sound  mathematical  foundation  as  well  as  developing 
methods  for  producing  such  vectors.  To  this  end,  the  analogy  with  contin- 
uous systems  has  been  invoked  extensively.  In  discussing  continuous  sys- 
tems, special  emphasis  has  been  placed  on  developing  criteria  for  the  selec- 
tion of  admissible  functions  capable  of  an  accurate  representation  of 
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the  motion  of  continuous  members,  as  well  as  on  developing  methods  for 
producing  such  admissible  functions.  In  regard  to  the  latter,  the  inte- 
gral formulation  of  the  eigenvalue  problem  (Ref.  5),  which  contains  the 
Green's  function  in  the  kernel,  proves  to  be  a useful  device  for  build- 
ing in  the  smoothness  and  boundary  conditions  requirements. 

Whereas  the  analogy  between  discrete  and  continuous  systems  makes 
physical  sense,  there  are  many  mathematical  aspects  requiring  close 
scrutiny.  For  the  most  part,  they  are  related  to  differences  between 
infinite  and  finite  vector  spaces.  In  some  ways,  the  formulation  for 
discrete  systems  is  simpler,  but  the  two  mathematical  models,  discrete 
and  distributed,  cannot  be  divorced  entirely.  Indeed,  one  must  consider 
the  origin  of  the  discrete  system  before  being  able  to  impose  proper 
requirements  on  the  smoothness  and  end  conditions  of  the  admissible  vec- 
tors. 

To  introduce  the  above  ideas,  a discrete  model  consisting  of  a 
rotating  lumped-parameter  cantilever  beam  has  been  used.  Such  a model 
is  of  special  interest,  as  it  can  be  regarded  as  simulating  a helicopter 
blade.  Because  the  continuous  counterpart  of  the  model  considered  is  a 
rotating  distributed-parameter  cantilever  beam  in  bending,  the  smooth- 
ness requirements  and  geometric  end  conditions  are  more  stringent  than 
those  for  a set  of  lumped  masses  simulating  a bar  in  axial  vibration. 

Once  again  drawing  on  the  analogy  with  continuous  systems,  ways  for 
producing  admissible  vectors  have  been  explored.  One  way  is  simply  to 
discretize  admissible  functions.  Based  on  the  recognition  that  the 
discrete  counterpart  of  the  Green's  function  in  the  integral  formulation 
is  the  flexibility  influence  coefficient,  another  way  is  to  use  the 
system  flexibility  matrix  to  generate  vectors  with  the  desired  smoothness 
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and  end  conditions. 

Preliminary  results  fully  justify  the  concept  of  admissible  vectors 
as  a valid  tool  in  dynamic  synthesis.  A comparison  of  eigenvalues  and 
eigenvectors  obtained  by  using  different  sets  of  admissible  vectors 
proves  quite  encouraging.  It  also  enhances  the  understanding  of  the 
problem  and  points  the  way  to  future  directions  of  the  research. 


where  X is  a parameter,  L is  a linear  homogeneous  differential  operator 
of  order  2p,  and  M is  merely  a function  of  the  spatial  variables.  Equa- 
tion (1)  must  be  satisfied  by  the  function  u at  every  point  of  an  open 
domain  D.  Associated  with  the  differential  equation  (1)  there  are  p 
boundary  conditions  that  u must  satisfy  at  every  point  of  the  boundary  S 
of  the  domain  D.  We  assume  that  the  boundary  conditions  are  of  the  type 

it 

BjU  = 0 , i = 1,2, . . . ,p  (2) 

where  are  linear  homogeneous  differential  operators  involving  de- 
rivatives normal  to  the  boundary  and  along  the  boundary  through  order 
2p-l.  The  eigenvalue  problem  consists  of  determining  the  values  of  the 
parameter  X for  which  there  are  nontrivial  functions  u satisfying  the 
differential  equation  and  the  boundary  conditions.  The  parameters  X are 
called  eigenvalues  and  the  corresponding  functions  u are  called  eigen- 
functions (Ref.  5). 

Quite  often  continuous  systems  lead  to  eigenvalue  problems  of  the 
type  defined  by  Eqs.  (1)  and  (2)  that  do  not  lend  themselves  to 
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closed-form  solution.  In  such  cases  we  must  be  content  with  approximate 
solutions  of  the  eigenvalue  problem.  We  shall  consider  the  Rayleigh- 
Ritz  method  for  obtaining  such  approximate  solutions.  Before  discussing 
the  Rayleigh-Ritz  method,  we  shall  present  certain  mathematical  prelim- 
inaries. 

Let  us  consider  a function  u that  is  continuous  in  the  closed  do- 
main D”D+S.  The  function  is  assumed  to  be  quadratically  integrable  in 

2 

the  Lebesgue  sense,  which  implies  that  the  integral  ||u||  = (u,u) 

= u“dD  exists,  where  ||u||  is  known  as  the  norm  of  u.  Our  interest 

lies  in  approximating  the  function  u by  some  other  function.  To  this 

end,  we  consider  the  sequence  of  functions  u^,  Uj,  •••,  u^,  ...  defined 

over  D.  Then,  if  lim  ||uu-u||  = 0,  it  is  said  that  the  sequence  Ujj 
N-*» 

converges  in  the  mean  to  u. 

Next,  let  us  assume  that  the  function  u can  be  represented  in  the 

OO 

domain  D by  the  series  u - £ a,i}>,,  which  is  convergent  in  the  mean, 

j“l  ^ ^ 

and  consider  the  partial  sum 
N 


“» ‘ £ Vj 


(3) 


Then,  the  set  of  functions  is  said  to  be  complete  if  it  is  possible 
to  find  an  integer  N and  a set  of  coefficients  a^,  a2»  ...»  a^  such  that 
for  any  £>0,  Hu-u^ll  < £.  The  question  remains  as  to  the  nature  of 


the  functions  <J)j  (j*l,2,...).  To  answer  this  question,  we  wish  to  delve 
a little  deeper  into  the  properties  of  the  differential  operator  L.  The 
operator  L is  said  to  be  self-adjoint  if  for  any  two  functions  u and  v 
in  the  field  of  definition  of  L 

(Lu,v)  * / vLudD  - / uLvdD  ■ (u,  Lv)  (4) 

D D 

Integrating  by  parts,  and  considering  the  boundary  conditions  of  the 
problem,  we  can  write 
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/ viudD  " / vApUdD  + / vAgUdS 


(5) 


D D S 

where  and  Ag  are  "two-sided"  symbolic  operators  symmetric  in  u and  v. 

We  should  point  out  that,  whereas  L is  an  actual  differential  operator, 

A is  an  operator  only  in  a symbolic  sense.  For  convenience,  we  shall 

use  the  notation 

[u,  v]|^  =*  J vAjjUdD  + J vAgUdS  (6) 

D S 

The  operators  L and  A are  said  to  be  positive  definite  if  the  inequality 
[u,u]A  > 0 is  satisfied. 

In  our  particular  case,  [u,u]A  defines  twice  the  potential  energy 
associated  with  a given  closed  domain  D.  Then,  the  quantity  |u| 
defined  by  I u | ■ [u,u],  is  known  as  the  energy  norm  of  u.  When  it  is 

A A 

clear  that  the  norm  involves  the  operator  A,  the  subscript  A can  be 
ignored.  The  sequence  of  functions  u^,  U2>  •••«  u^,  ...  converges  in 
energy  to  u if  lim  |cl^-u|  ■ 0.  If  Ujj  represents  the  partial  sum  (3), 

then  if  for  any  e>0  it  is  possible  to  find  an  integer  N and  a set  of 
coefficients  a^  a2>  ....  Sjj  such  that  lu-u^l  < e,  the  set  of  func- 
tions 0^ , <{>2>  •••»  <f>N»  •••  is  said  to  be  complete  in  energy. 

To  generalize  these  ideas,  let  L be  a positive  definite  operator 
defined  on  some  Hilbert  space  H.  But  the  scalar  product  (6)  also  de- 
fines a Hilbert  space.  The  Hilbert  space  defined  by  the  scalar  product 
(6)  is  called  complete  if  any  sequence  of  its  elements  satisfying  the 

condition  lim  * 0 has  a limit  which  is  in  the  space;  other- 

k,n-**> 

wise  it  is  called  incomplete.  In  general,  it  is  incomplete  and  must  be 
made  complete  by  adding  certain  new  elements  to  it.  The  new  Hilbert 
space  thus  constructed  is  called  the  energy  space  and  is  denoted  by  HA> 
Hence,  the  functions  must  belong  to  the  space  HA»  If  they  belong 
also  to  the  domain  of  existence  of  L,  then  [4)^,  ] * (M^,  ) - The 

question  is  as  to  what  distinguishes  the  functions  belonging  to  H from 
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those  belonging  to  H^.  It  is  clear  that  the  integration  by  parts  lead- 
ing to  Eq.  (5)  lowers  the  order  of  the  operator  involved  by  a factor  of 
two.  Hence,  if  L is  of  order  2p,  then  is  of  order  p,  so  that  the  re- 
quirements on  the  differentiability  of  the  functions  <J>^  are  lowered  ac- 
cordingly. More  importantly,  however,  the  energy  integrals,  Eq.  (6), 
take  automatically  into  account  the  natural  boundary  conditions,  so  that 
the  functions  need  satisfy  only  the  geometric  boundary  conditions. 
Such  functions  are  sometimes  referred  to  as  energy  functions.  In  more 
familiar  terminology,  we  shall  refer  to  functions  belonging  to  H as  com- 
parison functions  and  those  belonging  to  as  admissible  functions  (see 
Ref.  5). 


3.  Rayleigh's  Quotient  and  the  Rayleigh-Ritz  Method 

Let  us  assume  that  the  interest  lies  in  a self-adjoint  positive 
definite  continuous  system  whose  eigenvalue  problem  is  defined  by  Eqs. 
(1)  and  (2).  Let  v be  a comparison  function,  satisfying  all  the  boun- 


dary conditions  (2) , and  write 

2 LvlvdD 

u = R(v)  = ■ - ■ - 

/D *v  dD 


(7) 


where  R(v)  is  known  as  Rayleigh's  quotient  for  continuous  systems;  it  is 
always  positive  for  positive  definite  systems.  Then,  defining  a weighted 

z z 

norm  |]v||  ^ = ^V,V^M  = J Mv  dD  and  integrating  the  numerator  of  Eq.  (7) 
by  parts,  Rayleigh's  quotient  (7)  can  be  written  in  the  equivalent  form 

M2 


0)  - R(v) 


[v.v]^ 


(8) 


M INI  m 

The  advantage  of  Rayleigh's  quotient  in  the  form  of  Eq.  (8)  over  that  in 
the  form  of  Eq.  (7)  is  that  the  function  v in  Eq.  (8)  need  be  only 
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an  admissible  function.  In  either  form,  it  can  be  shown  (Ref.  5)  that 

Rayleigh’s  quotient  has  a stationary  value  in  the  neighborhood  of  any 

eigenfunction  u^r\  Moreover,  the  stationary  value  is  a minimum  in  the 

neighborhood  of  the  first  eigenfunction,  so  that  Rayleigh's  quotient 

2 

provides  an  upper  bound  for  the  first  eigenvalue  with  respect  to 

the  class  of  either  comparison  or  admissible  functions,  depending  on 

whether  Eq.  (7)  or  Eq.  (8)  is  used.  The  essence  of  the  Rayleigh-Ritz 

method  is  to  converge  to  the  desired  eigenvalues  by  constructing  a 

minimizing  sequence  to  represent  v. 

Let  us  consider  an  approximation  to  the  solution  of  the  eigenvalue 

problem  (1)  in  the  form  of  the  finite  series 
n 


l Vi 

i-1 


(9) 


where  <J>  (i-1, 2, . . . ,n)  are  known  functions  of  the  spatial  coordinates 

and  a^  are  unknown  coefficients  to  be  determined.  The  functions  (fK  are 
linearly  independent  over  the  domain  D and  are  either  comparison  func- 
tions or  admissible  functions  depending  on  whether  Rayleigh's  quotient 
is  in  the  form  (7)  or  (8) , respectively.  Following  a substitution  of 
Eq.  (9)  into  Rayleigh's  quotient,  the  coefficients  ai  are  chosen  so  as 
to  render  the  quotient  stationary.  It  can  be  shown  (Ref.  5)  that,  for 
Rayleigh's  quotient  to  be  stationary,  the  coefficients  a^  must  satisfy 
the  n algebraic  equations 


j 


Li 

I { C 4>  i * <f>j  1 a " Wai  3 0 


(10) 


where  the  form  (8)  of  Rayleigh's  quotient  has  been  used.  Note  that 
denotes  the  estimated  value  of  Rayleigh's  quotient,  in  which  n represents 
the  number  of  terms  in  series  (9).  Equations  (10)  represent  an  algebraic 
eigenvalue  problem  for  an  n-degree-of-freedom  discrete  system.  In  fact, 
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the  n-degree-of- freedom  system  can  be  regarded  as  the  result  of  imposing 
the  constraints 

an+l  = an+2  = •”  " 0 (11) 

on  the  inf initely-many-degree-of-f reedom  system.  Because  constraints 

have  a tendency  to  raise  the  stiffness  of  the  system,  the  estimated 

eigenvalues  A”  (the  square  of  the  natural  frequencies)  will  provide 

upper  bounds  for  the  true  eigenvalues  Xr> 

r =*  1,2, ...  ,n  (12) 


»n  . , 

A > a , 

r — r 


The  estimated  eigenfunctions  associated  with  the  estimated  eigenvalues 


are  obtained  by  introducing  the  coefficient  a^r^  in  Eq.  (9), 


» 


i 4r\ 

i=l 


(13) 


The  question  remains  as  to  how  the  estimated  eigenvalues  A^  and 


. (r) 


the  coefficients  a^  behave  in  the  limit.  Ee  adding  a term  to  the 


series  (9),  the  number  of  constraints  imposed  on  the  system  is  reduced 


by  one.  If  the  eigenvalues  Af  of  the  original  system  are  such  that 


i^<  A^<  ...  £ A“  and  if  the  eigenvalues  A“TX  of  the  system  obtained  by 


n+1 


reducing  the  number  of  constraints  by  one  are  such  that  _< 


.n+1 


^n+i»  then  by  the  inclusion  principle  (Ref.  3),  we  have 
»n+ 1 , .n  ^ .n+1  _ „ .n+1  ^ ,n  ^ .n+1  , , . . 

\ i\<*2  — * ’ * — An  — — An+l  (14) 

The  inclusion  principle  permits  us  to  conclude  that  the  estimated 
eigenvalues  tend  to  decrease  with  each  additional  degree  of  freedom.  At 
the  same  time  there  is  a new  eigenvalue  added  which  is  higher  than  any 
of  the  previous  ones.  The  estimated  eigenvalues  decrease  monotonically , 
as  the  number  of  degrees  of  freedom  increases.  Hence,  as  n-*»  the  esti- 


mated eigenvalues  A approach  the  true  eigenvalues  from  above. 


1 im  An  **  X 

r r 

n-wo 


1.2,... 


(15) 


■I 


8 


The  question  of  convergence  of  the  coefficients  a^  depends  on  the 

C r ) 

set  of  functions  <{k.  Determining  the  coefficients  a^  as  solutions  to 

the  eigenvalue  problem  (10)  is  equivalent  to  finding  the  best  approxi- 

(r) 


mation  to  the  true  eigenfunctions  u 
i.e. , 

(r)  (r) 


with  respect  to  the  energy  norm, 


I 


- u 


I A “ I’ 


^ £ a^r^<t>.|  =*  minimum 

i-1  1 1 


(16) 


Hence,  we  conclude  that  if  the  set  of  functions  <{>.  is  complete  in  the 


energy  space  H , u^  converges  in  energy  to  u^. 
a n 


Furthermore,  if  the 


set  of  functions  4^  is  minimal*  in  the  energy  space  HA,  then  limits  for 
n ( r) 

the  coefficients  a^  exist  and  have  the  values 


.<r) 


Urn  na^r) 

n-+oo 


1,2,. 


(17) 


Moreover,  if  the  system  of  functions  which  is  orthonormal  to  the  set  4k 
is  bounded  in  H^,  then  the  limiting  process  in  Eq.  (17)  is  uniform  with 
respect  to  i (Ref.  6). 

Implicit  in  the  above  discussion  is  the  need  to  approximate  a 
continuous  system  by  a discrete  one.  Indeed,  the  process  of  construc- 
ting the  series  (9)  is  called  discretization  because  it  reduces  the  pro- 
blem of  finding  an  infinite  set  of  continuous  functions  u to  that  of 
determining  a finite  number  of  constant  coefficients  a^  Of  course,  the 
object  is  to  obtain  an  accurate  approximation  to  the  lowest  eigenvalues 
and  eigenfunctions  while  taking  only  a small  number  of  n functions  in 
Eq.  (9).  The  question  remains  as  to  how  to  select  the  admissible 
functions  4>^  in  Eq.  (9). 


* A system  of  elements  of  the  Hilbert  space  H.  is  called  a minimal  sys- 
tem in  H,  if  the  deletion  of  any  one  of  the  elements  from  the  system  re- 
stricts the  span  of  this  new  set  to  a proper  subspace  of  the  space 
spanned  by  the  original  system  (Ref.  6). 
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4.  On  the  Selection  of  Admissible  Functions 

An  in-depth  discussion  on  the  rational  choice  of  the  functions 
can  be  found  in  Ref.  6.  The  discussion  in  Ref.  5 is  concerned  with  many 
aspects  of  the  problem  of  selection  of  the  functions  <p^,  including  the 
numerical  stability  of  the  approximate  solution.  However,  problems  of 
numerical  stability  are  of  no  concern  here,  because  quite  often  in  en- 
gineering it  is  possible  to  obtain  an  accurate  approximate  solution  of 
the  eigenvalue  problem  before  numerical  stability  becomes  a problem.  As 
a result,  we  will  be  concerned  only  with  obtaining  a set  of  suitable 
admissible  functions.  Hence,  ignoring  numerical  stability  aspects, 
criteria  for  the  selection  of  the  functions  0^  in  a closed  domain  D 
are: 

1)  Any  finite  set  of  admissible  functions  must  be  linearly  inde- 
pendent. 

2)  The  set  of  admissible  functions  must  be  complete  in  the  energy 
space  Ha  of  the  domain  D. 

In  addition,  as  mentioned  in  Sec.  3,  limits  for  the  coefficients  a^  in 
the  series  (9)  exist  if  the  set  of  admissible  functions  is  minimal. 

The  second  of  these  criteria  requires  that  we  be  able  to  identify 
a set  of  admissible  functions  that  is  complete  in  the  energy  space  H^. 

To  this  end,  the  following  statement  is  useful:  If  A and  B are  positive 
definite  operators  and  the  spaces  H^  and  contain  the  same  elements, 
then  any  set  that  is  complete  in  Rg  is  complete  in  HA-  Therefore,  by 
letting  B have  a simple  form  for  which  a complete  set  of  eigenfunctions 
is  easily  found,  we  can  use  that  set  as  a set  of  admissible  functions 
for  the  operator  A. 
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The  above  discussion  says  nothing  about  the  rate  of  convergence  of 
approximate  eigenvalues  and  eigenfunctions  for  a particular  set  of  func- 
tions satisfying  the  above  criteria.  In  fact,  the  approximate  eigen- 
solution  will  not  necessarily  converge  very  rapidly  to  the  exact  solu- 
tion. However,  in  practice,  we  are  concerned  with  finding  only  the 
lower  eigenvalues  and  eigenfunctions  and  relatively  simple  admissible 
functions  are  known  to  yield  rapid  convergence  to  the  lower  eigenvalues 
and  eigenfunctions  for  many  problems.  Therefore,  the  lack  of  analytic 
knowledge  about  the  rate  of  convergence  of  the  approximate  eigensolution 
for  a particular  set  of  admissible  functions  is  not  a deterrent  to  using 
the  Rayleigh-Ritz  method.  Indeed,  experience  has  shown  that,  by  compar- 
ing the  approximate  lower  eigenvalues  and  eigenfunctions  obtained  by 
using  several  different  sets  of  admissible  functions,  one  can  find  a set 
of  functions  which  yields  good  convergence. 

So  far  we  have  been  concerned  with  selecting  a set  of  admissible 
functions,  i.e.,  a set  of  functions  which  are  p times  differentiable  and 
satisfy  the  geometric  boundary  conditions . for  the  eigenvalue  problem 
(1) . We  shall  now  consider  a method  of  obtaining  not  just  admissible 
functions  but  comparison  function,  i.e.,  functions  which  are  2p  times 
differentiable  and  satisfy  both  the  geometric  and  natural  boundary  con- 
ditions of  the  eigenvalue  problem  (1).  The  method  is  based  on  the  inte- 
gral formulation  of  the  eigenvalue  problem  (1).  We  point  out  that, 
while  the  method  is  not  practical  for  finding  comparison  functions,  an 
analogous  method  for  the  discrete  problem  will  prove  very  useful. 

Indeed,  the  motivation  for  the  present  discussion  is  the  extension  to 
the  discrete  problem. 

The  eigenvalue  problem  (1)  can  be  written  in  the  equivalent  inte- 
gral form  (see  Ref.  5) 
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u(P)  - X /D  G(P,Q)M(Q)u(Q)dD(Q) 


(18) 


where  P denotes  the  position  of  the  point  at  which  the  displacement  u(P) 
is  measured  and  G(P,Q)  is  an  influence  function,  known  as  a Green's 
function.  The  Green’s  function  G(P,Q)  satisfies  all  the  boundary  condi- 
tions of  the  problem.  Moreover,  for  a self-adjoint  system,  the  Green's 
function  is  symmetric  in  P and  Q,  G(P,Q)  = G(Q,P).  Equation  (18)  can  be 
solved  by  iteration.  Denoting  by  Uq(P)  the  first  trial  function,  the 
first  iterated  function  u^(P)  is  given  by 

uL(P)  =*  JD  G(P,Q)M(Q)u0(Q)dD(Q)  (19) 

Then,  introducing  u^(Q)  in  the  integrand,  we  obtain  u^(P) , etc.  We 
are  not  concerned  here  with  iterating  to  the  actual  first  mode  but  with 
the  properties  of  the  functions  Uq(P)  and  u^(P).  Indeed,  if  Uq(P)  is 
any  function  which  is  continuous  in  the  domain  D,  then  u^(P)  will  be  2p 
times  differentiable  and  satisfy  all  the  boundary  conditions,  that  is 
u^(P)  will  be  a comparison  function.  In  addition,  for  a self-adjoint 
positive  definite  system,  if  Uq(P)  (i«l,2,...)  is  a set  of  linearly 
independent  continuous  functions,  then  the  set  of  functions  u^(P) 
(i-1,2,...)  will  also  be  linearly  independent. 

Let  us  consider  as  an  example  the  bending  of  a uniform  cantilevered 
beam  clamped  at  x*0.  In  this  case  Green's  function  G(P,Q)  reduces  to 
the  flexibility  influence  function  a (x,£)  where 

IT  2~  (x  " ^ - x 

a(x,£)  - (20) 

1 x^  1 

El  ~ (5  - 3 x)*  £>x 

in  which  El  is  the  bending  stiffness  of  the  beam.  The  eigenvalue  pro- 
blem has  the  form 


12 


1 


El  u""(x)  « Amu(x),  0 < x < L (21) 

4 4 

where  u""(x)  ■ d u(x)/dx  and  m is  the  mass  per  unit  length.  In  Eq. 
(21)  u(x)  satisfies  the  four  boundary  conditions 

u(0)  - u'(0)  - u" (L)  - u'"(L)  = 0 (22) 

Let  us  consider  the  two  linearly  independent  trial  functions 

Uq(x)  - C (23a) 

Uq(x)  - jr  (23b) 


Clearly,  these  functions  not  only  are  not  comparison  functions  but  they 
are  not  even  admissible  functions,  as  they  violate  the  geometric  boun- 
dary conditions  at  x*0.  Then, 

. L 

u.  (x)  * / a(x,Qm  u^(Qd£;,  i * 1,2  (24) 

x 0 u 

where 


uj(x) 


u1 2(x) 


CmL4  .xi  2x^1  A 
4EI  ^2  " 3 l3  6 l4; 


(25) 


mL 

6EI 


, *i  i *1 . __i 

L2  " 2 L3  20  L5 


) 


1 2 

and  it  is  easy  to  verify  that  u^(x)  and  u^(x)  satisfy  the  four  boundary 


conditions  (22)  and  are  linearly  Independent. 


5.  The  Discrete  Problem 

The  eigenvalue  problem  (10)  represents  a discrete  approximation  to 
the  continuous  eigenvalue  problem  defined  by  Eqs.  (1)  and  (2).  Quite 
often,  however,  engineering  problems  are  posed  directly  in  terms  of  an  n 
degree  of  freedom  discrete  system  yielding  the  eigenvalue  problem 

Ku  - AMu  (26) 

where  K and  M are  nxn  positive  definite  symmetric  stiffness  and  mass 
matrices,  respectively,  and  u is  an  n-dimensional  configuration  vector. 
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If  the  system  of  Eq.  (26)  represents  a mathematical  model  of  an  actual 
structure,  then  the  number  of  degrees  of  freedom  may  be  in  the  hundreds 
or  even  thousands.  Quite  often,  it  is  not  feasible  to  work  with  such 
large  systems  and  the  object  is  to  reduce  the  number  of  degrees  of  free- 
dom while  retaining  the  essential  dynamic  characteristics  of  the  original 
system.  The  approach  we  shall  consider  parallels  the  method  presented 
previously  for  continuous  systems.  But  first  we  shall  present  the 
mathematical'  concepts  for  discrete  systems  by  constructing  an  analogy 
with  the  concepts  presented  in  Secs.  2-4  for  continuous  systems. 

Let  H be  an  n-dimensional  vector  space  and  consider  the  n-dimen- 

2 


sional  vector  u in  H.  We  define  the  norm  ||u||  of  u by  ||u| 


(u,u) 


* u u.  If  we  wish  to  approximate  u by  some  sequence  of  vectors  u^,  U£, 

....  Ujj  (N  £ n)  in  H,  and  if  lim  || u^-u ||  = 0,  then  we  say  that  u^ 

N-^n 

converges  in  the  mean  to  u.  Let  us  assume  that  the  vector  u is  repre- 

~n 

sented  by  the  series  u =*  £ a . . which  is  convergent  in  the  mean,  and 

j”l 

consider  the  partial  sum 
N 


N < n 


(27) 


In  addition,  we  shall  define  the  energy  norm  for  discrete  systems  by 
2 T 

| u | * [u,u]  * u Ku.  The  sequence  of  vectors  u^,  u£,  ...,  u^  con- 

verges in  energy  to  u if  lim  larxl  " °*  If  Ujj  represents  the  partial 

N-m 

sum  (27)  then  the  linearly  independent  set  of  vectors  is  also  com- 
plete in  energy.  Of  course,  for  finite-dimensional  spaces  it  is  not 
hard  to  see  that  the  energy  space  HA  is  the  same  as  the  n-dimensional 

t 2 

vector  space  H.  If  we  now  define  the  weighted  norm  )|u||  M - (u,u)M  - 
T 

u Mu,  we  can  write  Rayleigh’s  quotient  for  the  discrete  system  (26)  as 
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0) 


R(v) 


[y,v] 

(^>M 


IyI 


(28) 


M 


where  v is  an  arbitrary  n-dimensional  vector.  The  Rayleigh  quotient 

(28)  has  a stationary  value  when  the  vector  v is  in  the  neighborhood  of 

(r) 

any  eigenvector  u (r»l,2, . . . ,n)  and  the  stationary  value  is  a minimum 
in  the  neighborhood  of  the  first  eigenvector. 

Next,  invoking  the  analogy  with  the  continuous  system,  we  select  a 
set  of  trial  vectors  4k  and  construct  a linear  combination 


~N 


l a <fr  (N  < n) 
i-1 


(29) 


where  the  4>^  are  known  vectors  and  the  a^  are  undetermined  coefficients. 
Substituting  Eq.  (29)  into  Eq.  (28),  we  choose  the  a^  so  as  to  render 
the  quotient  (28)  stationary,  yielding  the  N algebraic  equations 
N 


{[$t. 


(30) 


Equations  (30)  represent  an  N-dimensional  eigenvalue  problem,  so  that 

the  dimension  of  the  eigenvalue  problem  (26)  has  been  reduced  from  n to 

N 

N.  As  before,  the  estimated  eigenvalues  provide  upper  bounds  for  the 

N 

true  eigenvalues  X^,  A^  ^ Xf,  and  the  estimated  eigenvectors  are  obtained 
by  introducing  the  coefficients  a£r^  into  Eq.  (29) 

N 


* ai  $i 


(31) 


i-1 


N 


The  limit  properties  of  the  estimated  eigenvalues  Af  and  the  coeffic- 
ients a^r^  for  the  discrete  system  are  very  simple.  This  is  because  any 
set  of  n linearly  independent  vectors  4>^  is  both  minimal  in  the  energy 
space  and  complete  in  energy.  Hence,  if  the  vectors  4^  are  all  linearly 
independent,  then  mathematically  we  have 


15 


(32) 


lim  • 

N-hn  r 

11.  %(r) 

N-m  1 


V 


lim  | u 
N-m 


(t)--4r>i 


(r) 


.(r) 


and  the  limiting  process  for  the  coefficients  a/~'  is  automatically 
uniform,  i.e.,  all  the  limiting  values  are  reached  when  N=n. 


The  question  remains  as  to  how  to  choose  the  vectors  <p^. 


It  is 


obvious  that  for  discrete  systems  the  mathematical  analysis  does  not 
impose  restrictions  on  the  vectors  4*^  other  than  linear  independence. 

In  other  words,  the  analysis  does  not  guarantee  the  satisfaction  of  the 
boundary  conditions  and  smoothness  requirements  for  the  discrete  vectors 
4>^»  so  that  these  features  must  be  built  into  the  vectors  4>^  in  other 
ways.  In  fact,  if  the  discrete  system  represents  the  mathematical  model 
for  a continuous  structure  the  boundary  conditions  and  smoothness  require- 
ments are  already  taken  into  account  in  the  discrete  model.  This  sug- 
gests that  the  smoothness  properties  of  the  admissible  vectors  can  be 
obtained  from  the  discrete  model,  a subject  discussed  in  the  next  sec- 
tion. Our  interest  is  in  selecting  a small  set  of  N vectors  4^  which 
will  yield  accurate  results  for  the  lower  eigenvalues  and  eigenvectors 
of  the  n-dimensional  discrete  system. 


6.  The  Concept  of  Admissible  Vectors 

It  is  here  postulated  that  the  concept  of  admissible  functions, 
encountered  in  conjunction  with  the  mathematical  analysis  of  continuous 
systems,  can  be  extended  to  discrete  systems,  giving  rise  to  what  will 
be  referred  to  as  "admissible  vectors".  Essentially,  the  idea  is  that 
for  a particular  discrete  system  one  can  conceive  of  a simple  continuous 
system  that  is  similar  dynamically  to  the  discrete  system.  Such  a 
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continuous  system  admits  a set  of  admissible  functions  ({^(D).  By  choos- 
ing an  appropriate  set  of  n discrete  points  in  the  doi^ain  D of  the 
continuous  system,  we  can  construct  a set  of  n-dimensional  admissible 
vectors  <t>i  with  the  n entries  corresponding  to  the  values  of  the  func- 
tions ^(D)  at  the  n discrete  points.  This  process  guarantees  that  the 
admissible  vectors  are  in  some  way  "smooth".  Clearly,  if  the  admissible 
functions  are  linearly  independent,  then  the  admissible  vectors  con- 
structed by  "discretizing"  these  functions  will  also  be  linearly  inde- 
pendent. Moreover,  there  is  the  possibility  that  other  properties  of 
admissible  functions  carry  over  to  admissible  vectors.  The  validity  of 
the  concept  can  be  seen  by  considering  the  following  example. 

Consider  the  discrete  system  consisting  of  12  identical  masses 
connected  to  a support  and  to  each  other  by  12  identical  flexible 
massless  links  with  bending  stiffness  El  (see  Fig.  1).  The  support  is 
assumed  to  rotate  with  the  uniform  angular  velocity  ft.  We  shall  con- 
sider only  the  displacement  of  each  mass  from  the  plane  of  rotation. 
Taking  each  discrete  mass  to  be  equal  to  m,  the  12x12  mass  matrix  M is 
simply  the  12x12  identity  matrix  multiplied  by  m.  Calculation  of  the 
stiffness  matrix  K is  quite  a bit  more  involved.  In  fact,  it  will  prove 
convenient  not  to  calculate  K itself  but  to  calculate  its  inverse  K \ 
where  K * * A is  recognized  as  the  system  flexibility  matrix.  The 
entries  a^  of  the  matrix  A are  known  as  the  flexibility  influence  co- 
efficients, where  the  latter  are  the  discrete  counterparts  of  Green's 
function.  Obtaining  the  influence  coefficients  a^  is  no  simple  task  in 
the  case  at  hand.  Because  the  support  rotates  with  the  uniform  angular 
velocity  ft,  the  coefficient  a^  will  depend  on  all  the  coefficients 
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a.  (k-1,2, . . . ,12) . Hence,  for  each  j (j-1,2, ... ,12)  the  coefficients 
a.  (k-1, 2, . . . , 12)  can  be  shown  to  be  the  solutions  to  the  12  simultane- 

*^J 

ous  algebraic  equations 

n2.  3 i-1  . 4-1  12 

l {-|t(3i-l)  + 4[  l 2(i-r)  ] - 2(i-t)  ( J _ k)} 

R-  1 


r-1 


k-4+1 


2 3 i-1  12 

+ atj  U + ~r  ( 3 i—  1 ) + 2i  £ (i-r)  - 4 ( I k)]} 


n2,  3 
s)  h m 


2EI  3 

12 


r-1 


k-i+1 


k-i+1  J r=l 


i - 1,2,.. .,12 


(33) 


where  the  distance  between  each  mass  is  assumed  to  be  h.  Note  that  for 
each  j (j-1,2, ... ,12)  we  must  solve  a system  of  12  equations  in  the  12 
unknowns  a^  (k»l, 2, . . . ,12) . The  coefficients  a^  were  found  numeri- 
cally taking  h»L/12  and  - (l^)^  — and  are  exhibited  in  Table  1. 

m 

Hence,  we  are  concerned  with  the  discrete  eigenvalue  problem 

AmAu  - u (34) 

where  u is  a 12- dimensional  vector. 

The  eigenvalue  problem  (34)  can  be  solved  numerically  for  the  12 
actual  eigenvalues  A^  and  associated  eigenvectors  ur>  We  wish  to  com- 
pare the  actual  eigensolution  of  the  12th-order  eigenvalue  problem  (34) 
with  the  solution  of  the  reduced  eigenvalue  problem  (30),  which  in  terms 
of  the  flexibility  matrix  A takes  the  form 

l - $,V}a.  - 0 , i - 1,2 N (35) 

j-1  -1  1 

The  eigenvalue  problem  (35)  has  been  solved  five  times  using  different 
sets  of  vectors  The  first  solution  was  obtained  by  using  an  expan- 

sion in  terms  of  the  eigenvectors  of  the  nonrotating  discrete  system. 
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The  first  five  eigenvalues  and  eigenvectors  of  the  nonrotating  discrete 
system  are  displayed  in  Table  2.  The  eigenvalues  of  the  reduced-order 
rotating  system  using  a two-term,  three-term,  four-term  and  five-term 
expansiou  in  terms  of  the  nonrotating  eigenvectors  are  displayed  in 
Column  a of  Table  3.  Note  that  the  eigenvectors  of  the  nonrotating 
discrete  system  possess  all  the  smoothness  properties  of  the  discrete 
system.  A set  of  vectors  possessing  these  smoothness  properties  is  what 
we  will  call  a set  of  admissible  vectors.  Indeed,  the  bending  of  a 
cantilevered  beam  represents  a continuous  system  which  is  similar  to 
the  present  discrete  system  so  that  the  smoothness  properties  of  the  set 
of  admissible  vectors  correspond  to  the  differentiability  and  boundary 
conditions  requirements  imposed  on  the  set  of  admissible  functions  of 
this  similar  continuous  system.  The  admissible  functions  for  a canti- 
levered beam  are  functions  which  are  twice  differentiable  and  have  the 
displacement  and  slope  at  the  origin  equal  to  zero. 

The  second  solution  of  the  reduced  order  discrete  system  was  ob- 
tained by  using  an  expansion  in  terms  of  "discretized"  admissible  func- 
tions. We  take  as  a set  of  admissible  functions  of  the  cantilevered 
beam  the  functions 

<J>r(x)  - (x/L)r+1  r - 1,2,...  (36) 

and  construct  the  vectors  4>r  by  taking  as  the  ith  entry  in  the  12-dimen- 
sional vector  the  value  of  the  function  <p  at  y - ih  (i-1,2, . . . , 12) . 

r L 

Therefore, 

$r  - {(h) r+1,  (2h)r+1,  ....  (12h)r+1}T  (37) 

The  eigenvalues  of  the  reduced-order  system  using  the  vectors  of  Eq. 

(37)  are  displayed  in  Column  b of  Table  3 for  r-2,3,4,5,  respectively. 

Of  course,  it  is  easily  verified  that  the  functions  (36)  are  twice  dif- 
ferentiable and  that  <t>r (0)  - <f>'r(0)  » 0.  Hence,  the  vectors  (37) 
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possess  the  smoothness  of  the  functions  4>r* 

The  third  solution  of  the  reduced  discrete  system  was  obtained  by 

taking  discrete  values  of  the  functions 

4>r(x)  * 1 - cos(rrrx/L)  + y (-l)r+1(rTrx/L) ^ (38) 

at  the  points  — = ih  (i»l,2, . . . ,12)  as  the  entries  of  the  vectors  4>  . 

l r 

The  eigenvalues  are  displayed  in  Column  c of  Table  3 for  r»2,3,4,5,  re- 
spectively. Note  that  the  functions  (38)  are  not  only  admissible  func- 
tions but  comparison  functions,  i.e.  they  are  four  times  differentiable 
and  satisfy  all  the  boundary  conditions  4>r(0)  * 4^(0)  = <J>” CD  * 4>”'(L)  = 0. 

For  comparison  purposes,  the  fourth  solution  of  the  reduced  dis- 
crete system  was  obtained  by  using  a set  of  vectors  which  possess  none 
of  the  smoothness  characteristics  of  the  system.  The  motivation  is  the 
following:  if  we  use  a complete  set  of  twelve  arbitrary  linearly  inde- 
pendent vectors,  then  we  obtain  the  actual  solution  only  in  terms  of  a 
different  coordinate  system.  Hence,  the  question  arises  as  to  what  kind 
of  results  will  be  obtained  by  using  a reduced  set  of  these  vectors. 

The  following  vectors  have  been  used: 

- (1  1 1 1 1 1 1 1 1 1 1 1)T 

♦2  - (1  1 1 1 1 1 -1  -1  -1  -1  -1  -1}T 

<t>3  - {1  1 1 -1  -1  -1  -1  -1-111  1}T  (39) 

^•{111-1-1-1111-1  -1  -1}T 
♦5  - U l -l  -1-111  -l  -l  -l  l 1}T 
The  vectors  (39)  are  not  smooth  in  any  sense.  The  eigenvalues  of  the 
reduced  system  using  a two-,  three-,  four-,  and  five-term  expansion  in 
terms  of  the  vectors  (39)  are  displayed  in  Column  d of  Table  3. 

Next,  let  us  consider  the  discrete  analog  of  the  integral  formula- 
tion introduced  in  Sec.  4.  The  matrix  form  analogous  to  Eq.  (19)  is 
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uL  = K-1Mu0  = AMuq  (40) 

As  discussed  in  Ref.  7,  Eq.  (40)  is  the  basis  for  an  iterative  procedure 
to  obtain  the  actual  eigenvectors.  However,  because  our  interest  is 
merely  in  selecting  a set  of  vectors  which  will  yield  a good  representa- 
tion of  the  system  dynamic  characteristics,  it  is  not  necessary  for  us 
to  iterate  to  the  actual  eigenvectors.  Indeed,  we  shall  use  Eq.  (40) 
only  as  a way  of  imposing  the  smoothness  properties  of  the  discrete  sys- 
tem on  a set  of  arbitrary  linearly  independent  vectors.  If  the  system 
is  positive  definite  then  the  resulting  vectors  will  also  be  linearly 
independent.  Because  the  matrix  A has  built  into  it  the  system  smooth- 
ness characteristics  and  the  equivalent  of  boundary  conditions,  by  anal- 
ogy with  the  continuous  case,  the  vector  u^  will  possess  all  the  system 
smoothness  characteristics  regardless  of  the  form  of  Uq.  Hence,  the 
fifth  solution  of  the  reduced  discrete  system  was  obtained  by  using  the 
vectors  (39)  as  trial  vectors  on  the  right  side  of  Eq.  (40)  and  regard- 
ing the  resulting  vectors  as  admissible  vectors.  The  eigenvalues  for 
a two-,  three-,  four-,  and  five-term  expansion  are  displayed  in  Column 
e of  Table  3. 

The  actual  eigenvalues  of  the  12  by  12  system  are  displayed  in 
Column  f of  Table  3.  To  compare  the  actual  system  eigensolution  with 
the  eigensolutions  of  the  reduced  order  problem  (35)  for  the  five  dif- 
ferent 3ets  of  admissible  vectors  just  discussed,  the  error  in  the  ap- 
proximate eigenvalues  as  a percentage  of  the  actual  eigenvalues 
• N . 

M Aj  - X, 

ej  - ( - x )x  100  (41) 

has  been  calculated  for  two,  three,  four  and  five  term  expansions  for 


each  set  of  admissible  vectors.  The  results  are  displayed  in  Table  4, 
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where  Column  a of  Table  4 gives  the  percent  error  in  the  eigenvalues  of 


Column  a of  Table  3,  and  similar  statements  hold  for  Columns  b,  c,  d and 


e.  In  addition,  the  actual  system  eigenvectors  u as  well  as  the 


approximate  eigenvectors  u^  obtained  from  the  various  reduced  order 


eigenvalue  problems  have  been  calculated.  The  square  root  of  the  norms 


of  the  errors  of  the  approximate  eigenvectors 


*1-  HaS1*  - s<I)ll 1/2 


are  displayed  in  Table  5 for  two,  three,  four  and  five  term  expansions. 


where  Column  a of  Table  5 contains  the  norms  of  the  errors  of  the  eigen- 


vectors associated  with  the  eigenvalues  displayed  in  Column  a of  Table 


3 and  similar  correspondences  hold  for  Columns  b,  c,  d,  and  e.  Tables 


4 and  5 permit  us  to  make  several  remarks  about  the  relationship  between 


the  accuracy  of  the  approximate  eigensolutions  and  the  set  of  admissible 


vectors  used.  First,  Column  a of  Tables  4 and  5 reveals  that  in  using 


the  eigenvectors  of  the  nonrotating  discrete  system  to  reduce  the  order 


of  the  rotating  discrete  system  one  obtains  approximate  lower  eigenvalues 


and  eigenvectors  which  differ  only  slightly  from  the  actual  lower  eigen- 


values and  eigenvectors.  This  can  be  traced  to  the  similarity  of  the 


eigenvectors  of  the  nonrotating  discrete  system  with  the  eigenvectors  of 


the  rotating  discrete  system.  Hence,  the  choice  of  the  eigenvectors  of 


the  nonrotating  discrete  system  as  admissible  vectors  is  a good  choice 


provided  the  eigenvalue  problem  for  the  nonrotating  discrete  system  can 


be  easily  defined  and  solved.  However,  if  only  the  lower  eigenvalues 


and  eigenvectors  of  the  rotating  system  are  of  interest,  other  sets  of 


admissible  vectors  yield  equally  good  results  as  can  be  seen  by  inspec- 


tion of  Columns  b,  c and  e of  Tables  4 and  5.  Mote  that  using  the 


truncated  set  of  linearly  independent  vectors  (39)  as  admissible  vectors 


i 
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MMife*-, 


yields  very  poor  results  while  using  as  admissible  vectors  those  vectors 
obtained  by  imposing  the  system  smoothness  characteristics  on  the  vec- 
tors (39)  via  Eq.  (40)  yields  very  good  results.  This  can  be  seen  by 
inspection  of  Columns  d and  e of  Tables  4 and  5.  Moreover,  a comparison 
of  Columns  b and  e of  Tables  4 and  5 suggests  that  the  error  in  the 
approximate  eigensolution  obtained  by  using  "discretized"  admissible 
vectors  is  similar  to  the  error  in  the  approximate  eigensolution  obtained 
by  using  the  vectors  (39)  "smoothed"  by  Eq.  (40). 

The  degree  of  accuracy  of  the  approximate  eigenvectors  can  also  be 
seen  by  plotting  the  normalized  approximate  eigenvectors  u^.  Plots  of 
these  eigenvectors  for  two,  three,  four  and  five  term  expansions  are 
shown  in  Figures  2,  3,  4 and  5,  respectively.  The  approximate  eigenvec- 
tors are  obtained  by  using  admissible  vectors  as  in  Columns  a,  b,  c and 
e of  Tables  4 and  5 and  are  plotted  on  the  same  axes  as  the  actual 
eigenvectors.  It  is  apparent  from  Figs.  2-5  that  the  first  eigenvector 
is  always  very  accurate,  regardless  of  the  choice  and  the  number  of 
admissible  vectors  used.  Moreover,  as  we  increase  the  number  of  admissible 
vectors  used  in  the  expansion,  the  accuracy  of  the  second  and  third  ap- 
proximate eigenvectors  improves  regardless  of  the  choice  of  admissible 
vectors.  In  short.  Columns  a,  b,  c and  e of  Tables  4 and  5 along  with 
Figs.  2-5  permit  us  to  conclude  that  any  set  of  admissible  vectors,  i.e. 
vectors  possessing  the  system  smoothness  characteristics,  will  yield 
accurate  values  for  the  lower  eigenvalues  and  associated  eigenvectors. 

This  conclusion  should  have  far-reaching  implications  in  the  dynamic 
simulation  of  complex  structures  by  the  substructure  synthesis. 

7.  Summary  and  Conclusions 

In  this  investigation,  the  concept  of  "admissible  vectors"  for 
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the  simulation  of  discrete  structures  is  advanced  and  placed  on  a sound 


mathematical  foundation.  The  idea  is  to  simulate  the  motion  by  a sub- 
space of  smaller  dimnesion  than  the  dimension  of  the  actual  space,  i.e., 
to  truncate.  Admissible  functions  have  been  shown  to  be  sufficient  for 
the  simulation  of  distributed-parameter  members  in  a structure  consisting 
of  a given  number  of  substructures.  Admissible  vectors  are  intended  to 
play  the  same  role  for  discrete  substructures.  The  implication  is  that 
i set  of  vectors  satisfying  certain  smoothness  conditions  and  boundary 
conditions  will  yield  a satisfactory  representation  of  the  substructure 
motion.  Methods  for  generating  admissible  vectors  are  also  presented. 

To  introduce  and  validate  the  concept  of  admissible  vectors,  a discrete 
system  in  the  form  of  a lumped-parameter  rotating  helicopter  blade  has 
been  considered.  The  lower  eigenvalues  and  eigenvectors  obtained  by 
describing  the  system  in  terms  of  different  sets  of  admissible  vectors 
have  been  shown  to  converge  rapidly  to  the  actual  eigenvalues  and  eigen- 
vectors. These  preliminary  results  appear  most  promising. 

8.  References 

1.  Hurty,  W.  C. , "Dynamic  Analysis  of  Structural  Systems  Using  Compon- 
ent Modes",  AIAA  Journal,  Vol.  3,  No.  4,  1965,  pp.  678-685. 

2.  Meirovitch,  L.  "A  Stationarity  Principle  for  the  Eigenvalue  Problem 
for  Rotating  Structures",  AIAA  Journal,  Vol.  14,  No.  10,  1976,  pp. 
1387-1394. 

3.  Meirovitch,  L.  and  Hale,  A.  L. , "A  Rayleigh-Ritz  Approach  to  the 
Synthesis  of  Large  Gyroscopic  Structures",  Journal  de  Mecanique 
Appliquee,  Vol.  1,  No.  4,  1977,  pp.  333-356. 

4.  Meirovitch,  L.  and  Hale,  A.  L. , "Synthesis  and  Dynamic  Character- 
istics of  Large  Structures  with  Rotating  Substructures",  Presented 
at  the  IUTAM  Symposium  on  Dynamics  of  Multibody  Systems,  Munich, 
F.R.G.,  Aug.  29-Sept.  3,  1977. 

5.  Meirovitch,  L. , Analytical  Methods  in  Vibrations,  The  Macmillan 
Co. , New  York,  1967. 


24 


Mikhlin,  S.  G. , The  Numerical  Performance  of  Variational  Methods, 
Wolters-Noordhof f Publishing,  1971. 

Dong,  S.  B.  and  Wolf,  J.  A.,  Jr.,  "On  A Direct-Iterative  Eigensolu- 
tion  Technique",  International  Journal  for  Numerical  Methods  in 
Engineering,  Vol.  4,  1972,  pp.  155-161. 


25 


Table  1 - Matrix  of  Flexibility  Influence  Coefficients 


mm  ijg  MMl 


Table  2 

Eigenvalues  and  Eigenvectors  of  the  Nonrotating  Discrete  System 


A - diag  [.8776586  34.71710  273.9302  1058.188  2904.5644] 


6144262 *10~2 

. 3506042* 10-1 

. 8938545*10_1 

.1574684 

.2313173 

2363973*10-1 

.1198151 

.2648566 

.3836237 

.4309673 

5108426 *10-1 

.2244364 

.4026492 

.4044953 

.2050418 

8708854* 10-1 

.3215703 

.4199895 

.1600834 

-.2351863 

1302899 

.3883468 

.2947513 

-.1889837 

-.3938290 

1793706 

.4082550 

. 6761981 *10_1 

-.3947276 

-.9143334*10 

2330795 

.3724292 

-.1766866 

-.3102714 

.3186319 

2902570 

.2800194 

-.3431838 

-.8415515* 

lO-3  .3594727 

3498619 

.1374686 

-.3628411 

.3012934 

-.1365950*10 

4110011 

-.  4332252  •10-*1 

-.2171513 

.3626053 

-.3572553 

4729589 

-.2478191 

.5924247*10_1 

.1094921 

-.2424909 

5352285 

-.4624363 

.4007263 

-.3416881 

.2873619 

-1 
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Table  3 - Discrete  System  Eigenvalues 


No.  of 
Terms 

— 

Estimated 

Eigenvalues 

a 

b 

c 

d 

e 

f 

wmsm 

.3289033 

.3289754 

.3289252 

.3566144 

.3288961 

.3288935 

£ 

.4804341 

.5075637 

.4811161 

1.133973 

.4929387 

.4804187 

XjXlO-1 

.3288940 

.3288953 

.3288962 

.3559324 

.3288937 

.3288935 

3 

X2*10"2 

.4804338 

.4804287 

.4808920 

.5999323 

.4805295 

.4804187 

X3xlO-3 

.3099567 

.3498520 

.3111606 

.5434759 

.3166433 

.3099386 

A^xlO-1 

.3288936 

.3288935 

.3288940 

.3349127 

.3288935 

.3288935 

X2xlO-2 

.4804197 

.4804286 

.4804447 

.5578686 

.4804286 

.4804187 

A3x10-3 

.3099452 

.3105774 

.3111565 

.4531640 

.3109182 

.3099386 

A. xio-4 

.1129088 

.1361532 

.1139375 

.2147814 

.1286027 

.1129031 

Aj^xlO-1 

.3288935 

.3288935 

.3288936 

.3348560 

.3288935 

.3288935 

A2xlO-2 

.4804189 

.4804192 

.4804359 

.5514864 

.4804219 

.4804187 

5 

A3x10-3 

.3099391 

.3101618 

.3100193 

.4527830 

.3102866 

. 3099386 

A4xi0“4 

.1129044 

.1141827 

.1138691 

.1563050 

.1133845 

.1129031 

AjXIO-5 

.3019821 

.3855055 

.3066433 

.4084688 

.3156557 

.3019716 

a.  Nonrotating  eigenvectors  as  admissible  vectors 

b.  Admissible  vectors  as  per  Eq. 

c.  Admissible  vectors  as  per  Eq. 

d.  Admissible  vectors  as  per  Eq. 

e.  Admissible  vectors  as  per  Eq. 

f.  Actual  Eigenvalues 

J 


(37)  <*  = (jh)i+1) 

(38)  ($  - 1 - cos  jhiri  + y (-1) i+1(jhiri)2) 

(39) 

(AO) 
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Table  4 - Relative  Error  lr  the  Estimated  Eigenvalues  (in  7.) 


No.  of 
Terms 

Relative 

Error 

a 

b 

c 

d 

e 

9 

£1 

0.00296711 

0.02488588 

0.00962596 

8.42853092 

0.00077392 

e2 

0.00320669 

5.65029065 

0.14517522 

136.03847919 

2.60607359 

el 

0. 00014806 

0.00053885 

0.00081243 

8.22115504 

0.00004408 

3 

e2 

0.00314882 

0.00209765 

0.09851756 

24.87698487 

0.02307392 

£3 

0.00584250 

12.87784569 

0.39425222 

75.34951656 

2.16321827 

£1 

0.00001622 

0.00001260 

0.00016299 

1.83013831 

0.00000064 

A 

£2 

0.00020903 

0.00206151 

0.00542953 

16.12134874 

0.00207058 

£3 

0.00211595 

0.20608765 

0.39295927 

46.21087947 

0.31604395 

£4 

0.00503804 

20.59293695 

0.91621720 

90.23517802 

13.90534287 

£1 

0.00000258 

0.00000009 

0.00003613 

1.81288300 

0.00000017 

£2 

0.00004674 

0.00011852 

0.00358253 

14.79286722 

0.00068427 

5 

£3 

0. 00016512 

0.07201388 

0. 02603060 

46.08796571 

0.11228190 

£4 

0.00114977 

1.13339044 

0.85563648 

38.44174387 

0.42634579 

£5 

0.00349043 

27.66285680 

1.54706828 

35.26728795 

4.53158209 

Note:  Values  in  Columns  a,  b,  e correspond  to  admissible  vectors  as  in 

Table  3. 
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Table  5 - The  Norm  of  the  Error  in  the  Eigenvectors 


No.  of 
Terms 

Norm  of 
the  Error 

a 

b 

c 

d 

e 

61 

0.00547506 

0.01588587 

0.00986074 

0.28397162 

0.00279942 

6 

2 

0.00578902 

0.25104090 

0.03974361 

0.90269989 

0.17395070 

5 

1 

0.00121843 

0.00232567 

0.00285403 

0.27983097 

0.00066492 

3 

*2 

0.00572897 

0.00461749 

0.03204797 

0.47120213 

0.01554358 

*3 

0.00865562 

0.38862998 

0.06618765 

0.76453658 

0.17048115 

61 

0.00040298 

0.00035519 

0.00127729 

0.13463032 

0.00008023 

A 

62 

0.00145583 

0.00458223 

0.00740622 

0.38371777 

0.00459591 

63 

0.00484279 

0.04851585 

0.06602919 

0.61530317 

0.06034538 

54 

0.00871355 

0.49721320 

0.10629355 

0.94064578 

0.45201685 

61 

0.00016066 

0.00002935 

0.00060122 

0.13397904 

0.00004093 

52 

0.00068589 

0.00109361 

0.00600685 

0.36878886 

0.00262626 

5 

*3 

0. 00130988 

0. 02765621 

0.0 1633764 

0.61246579 

0.03441519 

*4 

0.00370327 

0.11857407 

0.10113619 

0.60093044 

0.07256947 

*5 

0.00778591 

0.58238777 

0.14596192 

0.68919375 

0.28400738 

Note:  Values  In  Columns  a,  b,  . . . , e correspond  to  admissible  vectors  as  in 
lible  3. 
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FIGURE  I * DISCRETE  SYSTEM  MODEL 
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b.  SECOND  EIGENVECTORS 


a,b,c,f 
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c.  THIRD  EIGENVECTORS 


FIGURE  4 * EIGENVECTORS  USING  A FOUR  -TERM  EXPANSION 
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FIGURE  5 - EIGENVECTORS  USING  A FIVE -TERM  EXPANSION 

(continued ) 


37 


